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Abstract 

This article studies the quasi-stationary behaviour of multidimen¬ 
sional birth and death processes, modeling the interaction between 
several species, absorbed when one of the coordinates hits 0. We study 
models where the absorption rate is not uniformly bounded, contrary 
to most of the previous works. To handle this natural situation, we 
develop original Lyapunov function arguments that might apply in 
other situations with unbounded killing rates. We obtain the expo¬ 
nential convergence in total variation of the conditional distributions 
to a unique stationary distribution, uniformly with respect to the ini¬ 
tial distribution. Our results cover general birth and death models 
with stronger intra-specific than inter-specific competition, and cases 
with neutral competition with explicit conditions on the dimension of 
the process. 
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1 Introduction 


This article is devoted to the study of quasi-stationary behavior of multi¬ 
type birth and death processes absorbed when one of the types goes extinct. 

More specifically, we consider a continuous-time Markov process {Xt, t > 
0 ) taking values in for some r > 1 , where we use the notations Z_|_ = 
{0,1, 2,...} and N = {1, 2, 3,..The transition rates of X are given by 


from n = (ni,..., Ur) to 


n + Cj with rate njbj{n), 


n — Cj with rate rij 


dj{n) + Y^l=iCjk{n 


for all 1 < j < r, with ej = (0 ,..., 0,1,0,... , 0 ), where the nonzero co¬ 
ordinate is the j-th one, 7 > 0, b{n) = {bi{n),... ,br{n)) and d{n) = 
(di(n),... ,dr{n)) are functions from to M!j_ and c(n) = (cjj(n))i<jj<r is 
a function from Z!j_ to the set of r X r matrices with nonnegative coefficients. 

In other words, the infinitesimal generator of the process X is defined 
for all bounded function / on Z!L and all n G Zj 


as 


Lf{n) = + ej) - f{n)]njbj{n) 




+ ej) - f{n)]nj 

j=i 


djin) + '^Cjkin)nk 


7-1 


\k=l 


( 1 . 1 ) 


This model represents a density-dependent population dynamics with r 
types of individuals (say r species), where bj{n) (resp. dj{n)) is the indi¬ 
vidual birth (resp. death) rate of an individual of type j in the population 
n, and Cij{n) represents the competition exerted by an individual of type 
j on an individual of type i in the population n. The global competition 
Ylk=i ejk{'n)nk felt by an individual of type i influences its death rate at a 
power 7 , which represents the strength of the competition. The larger 7 is, 
the stronger is the influence of the competition in large populations. The 
case 7 = 1 is very common in biology and is known as logistic (or competi¬ 
tive Lotka-Volterra) competition. Other values of 7 > 0 are also relevant in 
ecological applications. 

Note that the forms of the birth and death rates imply that, once a 
coordinate Xl of the process hits 0, it remains equal to 0. This corresponds 
to the extinction of the population of type j. Hence, the set 9 := Z!^ \ N*" is 
absorbing for the process X. Let us denote by tq its absorption time. Our 
goal is to study the process {Xt,t > 0 ) conditioned to non-aborption, and 
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in particular its quasi-stationary distribution, i.e. a probability measure a 
on N'" such that 

'^a{Xt G • I t < ra) = a, Vt > 0. 

More precisely, we shall give conditions ensuring the uniform exponential 
convergence of conditional distributions to the quasi-stationary distribution, 
independently of the initial condition. This means that there exist constants 
C, A > 0 such that 

\\F^{XtG-\t<Td)-a\\^y<Ce-^\ t>0, (1.2) 

where || • \\tv is the total variation norm and P(N^) is the set of probability 
measures on N^. This implies in particular the uniqueness of the quasi- 
stationary distribution. 

In particular, when (1.2) is satisfied, regardless of the initial condition, 
the quasi-stationary distribution describes the state of the population when 
it survives for a long time. One of the most notable features of quasi- 
stationary populations is the existence of a so-called mortality/extinction 
plateau: there exists Aq > 0 limit of the extinction rate of the population 
(see [16]). The constant — Aq is actually the largest non-trivial eigenvalue of 
the generator L and satisfies 

Pa(t < tq) = Vf > 0. 

Many properties can be deduced from (1.2), as the uniform convergence 
of < tq) to 7 ?(x), where rj is the eigenfunction of L corresponding 

to the eigenvalue Aq [6, Prop. 2.3], and the existence and the exponential 
ergodicity of the associated Q-process, defined as the process X conditionned 
to never be extinct (see [6, Thm. 3.1] for a precise dehnition). 

Quasi-stationary distributions for population processes have received 
much interest in the recent years (see the surveys [16, 22 ]). One of the most 
understood case concerns birth and death processes on absorbed at 0 : it 
has been shown in [ 20 ] that there are either zero, one or an inhnite contin¬ 
uum of quasi-stationary distributions for such processes. More recently, it 
has been shown in [15] that there exists a unique quasi-stationary distribu¬ 
tion for one dimensional birth and death processes if and only if ( 1 . 2 ) holds. 
This result has been extended in the recent paper [ 6 ] to birth and death pro¬ 
cesses with catastrophes (among other applications). The specific question 
of estimates on the speed of convergence to quasi-stationary distributions for 
one dimensional birth and death processes has been studied in [11, 8 ]. The 
multi-dimensional situation, which takes into account the existence of several 
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types of individuals in a population, is much less understood, except in the 
branching case of multi-type Galton-Watson processes (see [1, 17]) and for 
specific cooperative models with bounded absorption rate (see [ 6 ]). Another 
originality of the models studied in the present paper is the fact that the 
absorption rate is not uniformly bounded. To handle this natural situation, 
we develop original Lyapunov function arguments that might apply to other 
situations with unbounded killing rates. For results on the quasi-stationary 
behaviour of continuous-time and continuous state space models, we refer 
to [2, 14, 7] for the one dimensional case and to [18, 12, 3, 13, 10, 5] for the 
multi-dimensional situation. Infinite dimensional models have been studied 
in [9, 6]. Several papers studying the quasi-stationary behaviour of models 
applied to biology, chemistry, demography and finance are listed in [19]. 

We are going to prove (1.2) under two sets of assumptions. The first one 
(Section 2) considers stronger intra-specific competition than inter-specific 
competition (i.e. Cii{n) larger than Cij{n) for large |n| and for i ^ j). We 
make no particular assumption on the dimension r of the process and on 
the birth, death and competition functions. The second one (Section 3) 
considers a case of equal inter- and intra-specific competition (neutral com¬ 
petition), which leads to specific difficulties that we can solve for dimension 
r < 3 in the logistic case. The last section 4 is dedicated to a few extensions 
of our methods to other models. 

2 General birth and death functions with strong 
intra-specific competition 

The hrst case we study corresponds to the following assumptions, where |n| 
denotes ni for all n = (ni,..., Ur). 

Assumption 1. (HI) There exist constants h, d and c in (0, oo) and jdi > 
0 ,/32 G (— 00 , 1 ) with /3i -|- 7/^2 < 7 such that, for all n £ W and 
i G {1,.. . ,r}, 

0 < bi{n) < 6|n|^\ 0 < di{n) < J[n|^\ cu^n) > c|n|“^^. 

(H2) For all i £ {1,, r}, when jnj —>■ -|-oo, 

1 

Cn(n) > ^ Cjk{n) +— ^Cjjin), (2.1) 

l^j^k<r j=l 

where the notation f{n) ^ g{n) means that f{n)/g{n) —>■ -|-oo when 
\n\ -|-oo. 
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The first assumption is standard for population models, typically with 
/3i = 1^2 = 0. From the biological point of view, the second assumption 
corresponds to a stronger intra-specific competition than interspecihc com¬ 
petition. For example it holds if all the functions Cii{-) for 1 < f < r 
have the same order of magnitude when |n| ^ -|-oo and all the Cjj(-) 
for i ^ j are asymptotically negligible w.r.t. the Cjj('). The positivity of 
bi(n), Cii{n) ensures that the process is irreducible away from d (in the sense 
that Pn(-^i = m) > 0 for all n,m ^ d). 

Theorem 2.1. Under Assumptions 1, there exist constants C, A > 0 such 
that (1.2) holds true. 

We explain in Section 4 how to generalize this result to cases with killing 
or with multiple births. 

Remark 1. In fact. Assumption (H2) and the statement about Cii{n) in 
Assumption (HI) can be replaced by 

r / r \7 r / r \7 

E E ) Kj=i f ^ Cjk{n)nk I > 

j=i \fc=i / j=i \k=i / 

when |n| large enough, for some (explicit) constant Cr depending only on r 
and 7 . 

This allows to cover other biological settings than the strong intra- 
specihc competition. For example, assuming /3i = 0, it is satished when, 
for all i, Cij{n) = n\'^^ for at least some j ^ i for some (5 > 0 , and Cij{n) = 1 
otherwise. This corresponds to a situation of strong competition exerted 
on other species by large species (a kind of collective aggressivity against 
other species). Hybrid situations can also ht this assumption, for exam¬ 
ple when r = 2 , Ci 2 (n) = , C 22 {n) = nf and cii(n) = C 2 i(n) = 1 for 

some 5, 6' > 0. More complex structures of interaction can of course fit our 
assumptions. 

This also allows to cover classical biological cases, with comparable inter- 
and intra-specific competition, for example if 7 = 1 and Cij{n) = Cij inde¬ 
pendent of n, provided that the ca for 1 < f < r are not too small compared 
to the Cij j i 7^ j- 

This result and the one of Section 3 are consequences of the general 
criterion of [ 6 , Thm. 2.1], which gives a necessary and sufficient condition 
for (1.2) for general Markov processes. This condition is given by the two 
properties (Al) and (A 2 ) below: there exists a probability measure n on 
such that 
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(Al) there exist to,ci >0 such that for all x G N^, 

^x{Xto ^ -Ito < To) > cii/(-); 

(A2) there exists C 2 > 0 such that for all x (zW and t > 0, 

< Tq) > C2^x{t < Td)- 


Hence we only have to prove that Assumptions 1 implies (Al) and (A 2 ). 
Our proof of (Al) is based on Lyapunov functions and makes use of the 
following general inequality on conditional distributions. 


Proposition 2.2. Fix n G N*" and let /it(-) = lPn(Ait G • | t < tq). Let 
H —>• M_|_ such that LV is bounded from above on Then, for all 


t > 0, 


^it{V) - V{n) < f \^is{LV) - yis{V)iis{Llnr) 
Jo L 


ds, 


where the value of the integral in the r.h.s. is well-defined in (—oo, +oo] 


( 2 . 2 ) 

since 


fi 


fj.s{LV) - fisiV)fis{Llnr) 


ds < oo, 


where [x]_ = (— x) V 0 is the negative part o/x G M. 

Proof. Fix A: G N and let Tk '■= inf{t > 0 : \Xt\ > k}. We define as the 
process X stopped at time Tk, and denote by its infinitesimal generator, 
given by Lf f{n) = Lf(n)\\n\<k- Dynkin’s formula then entails 


E„F(Af) = V{n) + 


f 




L’^ViX^) 


ds. 


Letting k —)• +oo, Fatou’s lemma applied to both sides (mind that V is 
bounded from below while LV is bounded from above) imply that 

EnH(W) < V{n) + [ En [LV{Xs)] ds. (2.3) 

Jo 


Similarly, 


E„lp}’'(A'^^) — 1 + 


En 




ds. 


Jo 

Lebesgue’s dominated convergence theorem applies to the l.h.s. and the 
monotone convergence theorem to the r.h.s., which entails 


EnlN’'(-^t) 



[L 1 n.(X,)] 


ds. 
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Now, fix T > 0. We have 1 > Fn(i < tq) > ¥n{T < tq) > 0 for all t G [0, T], 
Note first that, since LV is bounded from above and V is non-negative, (2.3) 
implies that t ]Kn[LV(Xt)] G L^([0,T]). In addition, < 0 and 1^ > 0, 
so —E„(I/(W))E„(Lli!}r-(W)) > 0. Then, either 

[ E^(I^(X 5 ))E„(—Tlpjr (Xs)) ds =-|-oo, 

Jo 

and then ( 2 . 2 ) is trivial since utif) = E„(/(X 4 ))/E„(lN^(Xt)) > E„(/(Xt)) 
for all / > 0 , or 


f 


E„(I/(X,))E, 


(^ 5 )) ds ^ -j-oo. 


In this case, since nt{f) = E„(/(Xt))/E„(lN>-(-^t)) < 1E„(/(W))/P„(T < 
tq) for all / > 0 , we deduce from the fundamental theorem of calculus that, 
for all t G [0, T], 


f 


fJ,s{LV) - fIs(y)Hs{Lljqr) 


ds 


V{n) + /p*E„ [LV{Xs)\ ds 
^ /o [-^In^(^s)] ds 


Inequality (2.2) then follows from (2.3). 


□ 


Proof of Theorem 2.1. The proof of (Al) is based on the following bounded 
Lyapunov function: fix e G (0 ,7 — 7 / 32 ) and define for all n gW 


H .. 

= ^71+7’ 

i=i •' 

and Vs{n) = 0 for all n G 9 := Z!j_\N^. For all m,n gW such that \m\ < \n 
we have in particular the inequality 


i f -i- 

e V(l"^l + 1)^ 


\n\ + 1)' 


dx 


l\m\+l 




X 


/■l^l dx 

-1/ 

^ 1 

. 1 ^ 

l\m\ 

e \ 


\nY J 


. (2.4) 
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Step 1: Lyapunov function for the conditional distributions. 

Fix no € and let pt{-) '■= G ■ \ t < tq). It follows from (1.1) that 


p.t{LVe) = ht{n)Y 


+ lni=l 


nSN’’ i=l 

bi{n) 


nibi{n) 

1 (|n| + l)i+- 


m 


dii'^) + Cijin)nj 


(|n| + l)i+^ 


diin) + j y~]cij(n)nj 

i=i 


|n|i+= 


Vein) ^ (2.5) 


Fix n (zW and let i*(n) be (one of) the argmax of i i—)■ n* and let c*(n) := 
Ci*{n)i*{n)in). Then ni*(„) > |n|/r and 
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ni*{n) I I - 


c*(n)T'|nli+T' 


- 1+7 


Using Assumption (HI) and since diin) > 0, 


M,(iU) < X: M,(n) I - 1^,1. 


new 


Y = l 


2=1 


diin) + j ^Cij(n)nj 

1=1 


K(n)V. (2.6) 


In addition, 


r 

-HtiVe)l^tiLlw) < IlKiloo Y X] 

nCN’’ i=l 


diin) + 



Denoting a := 


c*(i,...,iF 




, the last two equations imply 


PtiLVe) - ptiVe)yitiLlnT) < Y Mn) 

new 


'^ + a— 


c*in)^\n\^-^ 

7 - 1+7 


+ 
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where, by (2.4) 




ni=l 


nSN’’ 


2 = 1 


di{n) + y~]cij(n)nj 


vi=i 


(||■^^.||oo-^.(n)) 




h|n|^i + (l]jViCii(n)|n| +Cii(n)^ 


nGN’" 

< ^ /rt(n) 
nGN’’ 


i=l 


e\ny 


'■d\n\^ 


1-^ C.^JnlT'-^ 

+ _Ii_LJ- 


1 

Cij (^) + -y X] 




2=1 


where C..y^r is a positive constant such that x\ + .. .+x 7 < C'.^^r( 3 ^i + .. .+Xr)'^ 
for all xi,..., Xr > 0 . 

Using Assumption (H2), we see that there exist constants B > 0 and 
C > 1 independent of n such that 


/rt(LK) - fit{Ve)MLlw) < Y, B{1 + \nf^-^) 

new ^ 


2 ri+T' / 


nSN’' ^ ^ 

\ riGN^ / 


where the last inequality follows from the fact that Pi < 'j — 7/^2 and e < 
7 - lp2- 

Step 2: Proof of (Al). 

Step 1 and Prop. 2.2 imply that 


< — 1 “ 

£ ./o 



^ X] j ds. 


This implies that, for any initial condition no, there exists s < l/(eC') such 
that 

\ny-"<^^-^Ps{n) < 2C^. 

new 

Let us dehne the set K = {n £ N'’, < 4C^}, which is hnite since 

e G ( 0,7 — 7 /^ 2 )• Using the previous inequality and Markov’s inequality, we 
obtain that Ps{K) > 1 / 2 . 
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The minimum p = m.\nx^K^x{Xu = Vu € [l/eC^2/eC)) is 

positive because K is finite, X is irreducible away from d and the jumping 
rate from (1,..., 1) is finite. Hence 

/i,(p.(X„ = 1 ), Vug [l/eC, 2 /eC))) > | > 0 , 

Using the Markov property, we deduce that 

/^ 2 /(£C){(l! ■ ■ ■ , 1 )} > Ais(lP-(-^ 2 /(eC)-s = ( 1 ) • • • ! 1 )) 

>/i,(P.(X„ = 1 ), Vug [l/(eC),2/eC)) 

>f>0. 

Since p does not depend on the initial distribution of the process, we deduce 
that (Al) is satished with v = to = 2/(eC') and ci = p/ 2 . 

Step 3: Proof of (A 2). 

The same calculation as in (2.6) shows that, for all n G N^, 

rr.H < . 

Now, using Assumption (HI), we have c*(n)'^|np“^ > with 

7 —£— 7/32 > (/3i—£)V0. Hence, there exist two positive constants Ci, C 2 > 0 
such that 

LVe{n) <Ci- , Vn G NT 

Since V/ is bounded and LV/ is bounded from above, we deduce from 
Dynkin’s formula as in the proof of Proposition 2.2 that, for all A: > 1, 

< K(n) + E ^ LVeiXt)dtj 

< i + (Ci - E„ {T{\m\<k} A ra) , 

where is the hrst hitting time by Xt of the set {m G : \m\ < k} 

Since V£{Xr^^^^^f,yATg) is almost surely non-negative, we obtain 

sup E„ (Ts\m\<k} A Ta) - 0 . 

fc-i>oo 

Using the same argument as in [ 6 , Eq. (4.6)], we deduce that, for all A > 0, 
there exists A: > 1 such that 

sup E„(e^Ul-l<''>^"'a) < + 00 . (2.7) 
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Let us now denote by A the total jumping rate from (1,... , 1). We choose k 
such that (2.7) holds for this constant A. Defining the finite set G = {m € 
I \m\ < A:}, we thus have 

A := sup < oo. (2.8) 

nGN’’ 

The irreducibility of X and the finiteness of G entail the existence of a 
constant C > 0 such that 

supPn(t < Tg) < G inf P„(t < rg), Vt > 0. (2.9) 

n€G 

For all n € N'", we deduce from Chebyshev’s inequality and (2.8) that 

Pn(A < To A Tg) < Ae~^*. (2-10) 

Using the last two inequalities and the strong Markov property, we have 


Pn(A < Tg) = Fn(t < tg A Tg) + Pn(TG A Tg < t < Tg) 

< Ae~^^ + / sup Pm(t - s < Tg)Fn{TG A Tg ^ ds) 

Jo m&G\J{d} 

< Ae~^* + C [ Pxo(A - s < Ta)P„(TG ATg € ds), 

Jo 

where xq := (1,..., 1). Now, by dehnition of A and by the Markov property, 

Pxo(A - s < Tg)e~^^ = P3;o(t - s < Tg)Fa;o{Xu = Xo,Vu € [0, s]) 

< Pxo(t < Tg). 


Hence 


lPn(A < Tg) < Ae 


+ CFxq 



e^'^'^niTG A Tg € ds). 


We hnally deduce that 

^n{t < Tg) < AF^^{t < Tg) + GAF^^{t < Tg). 

This entails (A2) for v = 5 xq = 


□ 
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3 Birth and death processes with neutral compe¬ 
tition in not too large dimension 

The second case we study corresponds to the following assumptions. 

Assumption 2. (H3) There exist eonstants h, d in (0, oo) and /3i G [ 0 , 7 ) 
such that, for all n G and i G {1,... , r}, 

0 < bi{n) < 6|n|^^, 0 < di{n) < 

(H4) There exists a constant c > 0 such that, for all i,j G {l,...,r}, 
Cij{n) = c (^neutral competition^. 

Theorem 3.1. Under Assumptions 2 and ifr < 1 + 67 , there exist constants 
C, A > 0 such that ( 1 . 2 ) holds true. 

Remark 2. In the classical logistic case 7 = 1 , this gives the existence and 
exponential convergence in total variation to the quasi-stationary distribu¬ 
tion up to dimension 3. The larger 7 is, the larger the dimension r can be 
taken. 

Proof of Conditions (Al) and (A2) under Assumptions 2. We proceed as in 
the proof of Theorem 2.1. Fix e G ( 0 , 7 ). First, we deduce from (2.5) and 
from the fact that > n* — 1 that 



In addition, 


r 


-p.tiVe)pt{Llw) = Pt{Ve) ^ Pt{n) ^ lni=i [di{n) + c^\n\^] . 


neN’' i=l 
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Hence 


Ht{LVe)-nt{Ve)fJ.t{Llnr) < ^ fit{n) 
where 


nCN’' 


b\nf^ 


rc' 


n 


l+e—7 


+ ^, 

(3.1) 


H := ^ fitin) ^ l„,=i [di{n) + c'^|nP] (//t(14) - 14 (n)) 

nGN’’ 


i=l 


< Mn) 

new 


dr\nf^{\\Ve\\oo - Ve{n)) + c^r 


7 .^ 1+7 


+c'^|nP(r - 1) intiVe) - Vein))+] , 


(3.2) 


where the last two terms are obtained by distinguishing between the cases 
where n = (!,...,!) (and hence |n| = r) and n ^ (!,...,!) (and hence 
I]lni=i < r - 1). Now, defining We{n) := ||14||oo - 14(n), it follows 
from (2.4) that 

Y d't{n)\n[^ if^tiVs) - K(n))+ 


neN’' 


< 


- 7 /e 


/ii(n)lH,(n)-^/^ (Wein) - MWe)). 


nSN’’ 


-^Y f^tin)We{n)^ 


tiGN’’ 


MWs) 

We{n) 


Defining 


we obtain 


/ \ h't(ll^e) w TNjr 

We(n) := ^ ^ , Vn € N , 


Y d't{n)\n[< intiVe) - 14(n))+ < 
neN’' 

Now, it is elementary to check that 


^ ^ ^ / 7-£ 
X £ (1 — X) < — - 

7 V 7 

and Holder’s inequality implies that 


We{n) 

Y.new Mn)we{ny^''~^{l - We{n))+ 


1 < ( We 


eT'/" 


, Vx > 0 


e /7 
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Hence, it follows from (2.4) that 


Y1 (MVe) - Vein)). 


< 


/ii W, 


nCN’’ 


Since 


je 


7/e-l 


7 -e 


<- 1 -- 

7 V 7 


<■-1 


(ki + 


neN’" 


lim (1 — x)^ ^ . 

^ ^ e' 


under the assumption that r < 1 + 07 , we can hnd e > 0 small enough such 
that, for some <5 > 0 , 

- 1 ) (/^t(K) - Ve{n))^ < (1 - 5) J]] (|n| + 

nEN’’ riEN^ 

Combining this with (3.1) and (3.2), since /3i — e < 7 — e and e < 7 , there 
exists a constant C > 0 such that 


Ht{LVe) - fJ.t{Ve)fJ.t{Llw) < ^ lit(n) fb + —"j |n|^i ^ + c 


,7|„|7-£-l 


nEN' 

"y l“h'y 

+^2^^-c^'lnp-^ + c^{l - 6){\n\ + 1)^'“ 


e 


c^6 


' Mn)- 


nSN’’ 


The rest of the proof is the same as for Thm. 2.1. 


□ 


4 A few extensions to other models 

The method that we used is based on a Lyapunov type argument to prove 
conditions (Al) and (A2). This method is general enough to apply to a wide 
range of other models. We give here a few simple examples for which the 
exponential convergence of conditional distributions can be proved following 
the same arguments. 

One dimensional birth and death processes with catastrophes. We 

consider a standard birth an death process on Z_|_ with birth (resp. death) 
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rate bn (resp. dn) from state n G Z_|_, with do = bo = 0 and bn,dn > 0 
for all n G N. This process is absorbed at 9 = 0. Moreover, we add a 
catastrophe rate a„ > 0 of jump from any state n G N to the absorption 
point d. Such models have been studied in [21, 6] with an assumption of 
uniformly bounded catastrophe rate, which we relax here. 

In this example, we restrict for simplicity to the logistic cases, where 
there exist constants 6 > 0, c > 0 and 6 G (0,1) such that, for n large 
enough. 


bn < nb, dn > n^c and an < den. (4.1) 

This simple situation allows a similar computation as in Section 2 (with 
7 = 1 and /3i = /32 = 0). Of course, the arguments can be easily adapted to 
any other cases where explicit Lyapunov functions are known. 

Proposition 4.1. Assume that (4.1) holds for n sufficiently large, then 
there exist constants C, A > 0 such that (1.2) holds true. 

Proof. We use the same Lyapunov function Vg as in Section 2 and we obtain 
MLV,) = Y, l^tin) ( . 

neN VV ^ ; / 

< Y ~ - anVe{n)) . 

neN 


In addition 


-fItiVe)nt{LlN) < IIKII oo E Utin) (ln=lCl +an) ■ 

neN 

Hence there exists a constant C > 0 such that 

Ht{LVe) - Ht{Ve)nt{Lln) < C - //t(n)n^“^ + ^//t(n)a„ (||I4||oo - Ve{n)) 

7T.EN tiEN 

< C - cY ^1 - - V 

Taking e G (5,1), we can conclude as in Section 2. □ 


Multi-dimensional birth and death processes with catastrophes. 

We now study the multi-dimensional case with catastrophes, for which, as 
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far as we know, no result on quasi-stationary distributions is known. We con¬ 
sider the same model as in Section 2 with an additional jump rate a{n) > 0 
from any state n E N*" to d. The next result can be proved by an easy com¬ 
bination of the arguments of the proofs of Theorem 2.1 and Proposition 4.1. 

Proposition 4.2. Under Assumption 1 and the assumption 

a{n) < Cii{n)\n[<, 

there exist constants C, A > 0 such that (1.2) holds true. 

Multi-dimensional birth and death processes with multiple births. 

We consider the same model as in Section 2 with an additional feature: we 
allow multiple progeny at each birth time. To do so we consider, for all 
n E N*”, a probability measure 

Pn ■— ^ ^ Pn,khk' 
fceZ!;. 

Then, when a birth occurs (at rate b{n)) in a population n, the new state 
of the population is n -|- A: with probability Pn,k- A one-dimensional case has 
already been studied in [4]. 

Proposition 4.3. Under Assumption 1 and the assumption 
M := sup ^ \k\Pn,k < oo, 

there exist constants C,X> 0 such that (1.2) holds true. 

Proof. The only term which differs from the proof of Theorem 2.1 is the 
birth term in given by 


^ ^ Pn,k 7T+F 

nGN’' i=l j=\n\+l 

\k\ 

< ^ ^ -r^7j^Pn,k 

neN’" i=l ke.'L']^ ' ' 

< M ^ pLt{n)b\n\^^~^. 

Thus, we obtain a similar bound as in Section 2 for this term, and the proof 
can be completed as there. □ 
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